Constraints on a f{R) gravity dark energy model with early scaling evolution 
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The modified gravity witii f{R) = R^'^^ (e > 0) allows a scaling solution where the density of 
gravity sector follows the density of the dominant fluid. We present initial conditions of background 
and perturbation variables during the scaling evolution regime in the modified gravity. As a possible 
dark energy model we consider a gravity with a form f{R) — R^^"^ + qR~" ( — 1 < n < 0) where 
the second term drives the late-time acceleration. We show that our f{R) gravity parameters are 
very sensitive to the baryon perturbation growth and baryon density power spectrum, and present 
observational constraints on the model parameters. Our analysis suggests that only the parameter 
space extremely close to the ACDM model is allowed. 

PACS numbers: 98.80.-k, 95.36.+X 



I. INTRODUCTION 

The current accelerated expansion of the universe has 
been theoretically explained by either the cosmological 
constant or some kind of dynamical dark energy in the 
context of field or modified gravity (see Ref. [l[ for recent 
reviews). 

In the modified gravity such as f{R) gravity 0] (see 
Rcfs. Q for reviews and 0-13 for recent works), back- 
ground evolution sometimes confronts with a numeri- 
cal difficulty. The background evolution is not easily 
handled numerically during the early radiation domi- 
nated era in the modified gravity with a functional form 
f{R) = i? + /de(^), where /de drives the late-time ac- 
celeration. The cosmologically viable forms of /de pro- 



posed so far are qR ' 
Sawicki Q, Starobinsky 



and forms given by Hu & 
and so on [lol . [Tl| . In all 

models, the second term /de becomes extremely subdom- 
inant compared with R in the early radiation dominated 
era so that the f{R) gravity effectively goes over into 
the Einstein gravity; in general, however, the evolution 
could be more complicated, see [l^. Since the quantity 
F = df /dR becomes extremely close to unity in such a 
situation, evolving a differential equation like Eq. ([T} or 
(21) below is sometimes not numerically feasible. 

Adopting the modified form f{R) = R^~^^ + /de with 
small positive e together with appropriate initial condi- 
tions we can evade this numerical problem (here we use 
the Planck unit with SttC? = 1 = c). It is known that 
the first term R^~^'^ which is dominant in the early epoch 
allows the density of gravity sector to follow that of dom- 
inant fluid (scaling evolution) [11 . 

We are motivated to study the case in order to inves- 
tigate the observationally allowed regions with qualita- 
tively different evolution available in our case of f{R) 
gravity. By considering R^^'^ term, however, the grav- 
ity with /de = — 2A does not go over into the Einstein 
gravity in recent era. We will still consider values of e 
which is likely to be excluded by the solar-system test 
because with vanishingly small e the system of equations 



cannot be handled numerically due to limited numerical 
precision in the early era. Though, wc will show that for 
smallest value of e we considered, the cosmological evo- 
lution we study is numerically similar to the evolution in 
Einstein's gravity with cosmological constant. 

In this paper, wc present initial conditions of back- 
ground and perturbed variables during the scaling regime 
in this gravity. Using the initial conditions for scal- 
ing evolution we present background evolution, mat- 
ter (density) and cosmic microwave background (CMB) 
anisotropy power spectra, and perturbation growth in the 
gravity with /de = qR~^ ■ We show that the CMB power 
spectrum is not sensitive to the model parameters, and 
explore the viable parameter space constrained by the 
type la supernova (SNIa), matter power spectrum, and 
the future perturbation growth factor observation. 

Throughout this paper wc assume spatial flatness 
{K = 0). Notations and the basic set of equations in 
f{R) gravity are summarized in Ref. p^ . 



II. BACKGROUND AND PERTURBATION 
EQUATIONS 

The background evolution in the f{R) gravity is de- 
scribed by (see Eqs. (43), (57) and (59) of Ref. [ll) 

F + 3HF + i(2/ - FR) = i(/i,„ ~ 3p,„), (1) 

where H = d/a, a{t) is the cosmic scale factor, a dot 
represents a derivative with respect to the cosmic time 
t, and firn and Pm are collective density and pressure 
for radiation (i?) and matter (M): fim — fJ-R + fJ-M and 
likewise for pm- Using F = F^nR with F^n = dF/dR, 
Eq. ([1]) is transformed into a differential equation for R 
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3HR- 
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3F, 



■{(Im-SPm), (2) 



which is what wc actually solved numerically. To evolve 
this equation wc need H{t). For spatially flat Robertson- 
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Walker metric, we have (sec Eqs. (43) and (57) of Rcf. 

El) 

H^ = J^ (m™ + Fmx) , = 1{FR- /) - 

(3) 

where fix indicates the energy density of f{R) gravity 
sector (hereafter X-coniponent), and R = 6{2H^ + H). 
With the definition of conventional density parameters, 
= /^i/(3-ff^) {i — R, M, X), we have a relation (17 + 
flM)/F + fix — 1 from Eq. ([3]). The equation of state 
of the X-component is defined as wx = Px / l^x , where 



px is the pressure of the X-component given by (see Eq. 
(57) of Rcf. [H) 

Fpx = -^iFR-f) + F + 2HF. (4) 

The perturbation equations in the CDM-comoving 
gauge (CCG), which sets perturbed velocity of cold dark 
matter is set to zero {vc = 0) as the temporal gauge (hy- 
persurfacc) condition, are presented in Eqs. (86)-(88), 
(66), and (67) of Ref . Using SR as the perturbation 
variable, we have 
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3(5pm) , (5) 



(6) 



C = (l + -) (|-^-^ 
^4 + (1 - 3w) Uu, = 



1 + w aH 



(7) 
(8) 
(9) 



where Sfim and dpm are collective perturbed density, 
pressure, respectively. In our full numerical treatment 
the collective fluids includes the CDM, baryon, radiation 
(photons and neutrinos), etc; the radiation parts include 
photons handled by Boltzmann equations or tight cou- 
pling approximation. In Eqs. ([5]) and 0, as an exam- 
ple, we present a fluid with a constant equation of state 
parameter w = Pw/ fJ-w] — Sfiw/ fi^ is the density con- 
trast and is the perturbed velocity of the dominant 
fluid. Here a prime indicates a derivative with respect 
to In a {F' = dF/dlna), k is the comoving wave num- 
ber, and K is the perturbed expansion of normal frame 
vector. As we consider the CCG we need to include the 
CDM component even in the case it is subdominant. In 
case it is subdominant, we have at least a sixth order 
differential equation as presented in Eqs. ([5]) -([HI). 



III. INITIAL CONDITIONS FOR SCALING 
EVOLUTION 

In the early era, let us consider a functional form 

/(i?)=«i?i+% (10) 

where e is small positive constant; for generality we in- 
troduced a coefficient a which is unity in our unit. The 
gravity of pure power-law form allows scaling evolution 
in which the density of X-component follows that of 
the dominant fluid (Bj, and the corresponding poten- 
tial in the Einstein frame is a pure exponential poten- 
tial |13j . Here we put an ansatz that Ffix evolves as 
the dominant ideal fluid with constant equation of state 
{w ^Pw/y^w = Spw/Sp-w) as 



Ffix = \{FR-f)-3HF = Afi^, (11) 

where A is the constant density fraction of X-component 
relative to the dominant fluid. By combining Eqs. ([Ij, 
([3]), and (|TT|) . we can derive 



3FR 

F_ _ 
~F ^ 



-/ 
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= 2Ai„[(l 
3u;)[(l4 



-iw) 
A)f- 



M2- 

FR] 



f ^(H 



3w)FR 



(12) 

(l + ^)(/-3Fi?) ■ ^^^^ 

Equations and (|13p provide the scaling initial condi- 
tions for F and F' , respectively. Noting that 3FR — / cx 
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-3{l+w) 



in Eq. (|12p and specifying the form of /(i?) in variables, 



Eq. ((TU]), wc obtain an exact expression for R 



R — Hf, 



6r2^o[(l-3w) + (2-3w)A] 



a(2 + 3e)(a/ao)3(i+'") 
and its time derivative 



l + e 



1 



(14) 



(15) 



where fl^o = ii^o/{3Hq) is the density parameter of w- 
fluid at the present epoch (indicated by the subscript 0). 
Inserting Eq. (|15p into Eq. ()13|) . we determine 



e(7 + 10e) + 3e(l + 2e)«; 
2-3e-8e2-3e(l + 2e)w' 



(16) 



Equations p^ -pB )) can be used as initial conditions 
for the scaling density evolution of X-component in the 
f{R) = i?^"*"^ gravity. The scaling behavior is possible for 
general constant w value of the dominant fluid. As we set 
the scaling initial condition, the density of X-component 
follows (scales) the dominant fluid component even for 
changing w value of the dominant component; for exam- 
ple, from radiation dominated era to the matter domi- 
nated era. 

Now, for the perturbed initial conditions, applying the 
same method used in Ref. [l6| that presents initial condi- 
tions of perturbed variables during the scaling regime of 
the cosmology based on a minimally coupled scalar field, 
we find initial conditions of perturbed variables in our 
f{R) gravity. In the large-scale limit (^ 0), with the 
help of Eqs. (O, ([TU)) -([TC 1) . all the background-related co- 
efficients in Eqs. ([S])-® are expressed in terms of e and 
w alone during the scaling regime. We have solutions for 
Eqs. dl])-® in the form 



(17) 



where 



l-2e + 3w 




3 fl + w 
2 I l + e 



3e + 3e(l + 2e)u; ± 



16e 



(18) 



4e(l + e) 

- 316^ + 160e^ + 2566* + 3e(16 - 22e 

- 28e^ + 6Ae^)w - 96^(7 + 12e - Ae^)w^ 



The solutions have been obtained by MAPLE software of 
version 11 17[. By choosing the first one as the growing 



l~2e + 3w 
HJ " (l + e)(H-u;)' 
1 - 4e - 3{1 + 2e)w 



6R _ 

~R ~ 1 + 7e - 12e2 - 3(1 - e)(l - 2e)w^ 
5R\' l-2e + 3w f 5R^ 



(19) 



1 



R 



where C is the initial amplitude. 

Solutions similar to Eq. ^TE\i were presented in Ref. 
[l^ . Compared with the solutions in Ref. [3] which 
were based on the comoving gauge of the dominant fluid 
(the uj-fluid), our solutions are based on the CCG gauge 
condition which is the CDM comoving gauge. A compar- 
ison of the two results shows that the growing solution 
coincides and the other three solutions differ for w ^ 0. 
We are interested only in the growing solution, and the 
solutions in Eq. (fT9)) coincide with the one in Ref. (l^ . 



IV. 



A DARK ENERGY MODEL AND 
CONSTRAINTS 



In this section, we present numerical evolution of back- 
ground and perturbed quantities in a dark energy model 
based on the f{R) gravity with early scaling era. We 
consider a gravity with a form 



f{R) = R 



l + e 



qR- 



(20) 



mode solution, we get initial conditions of the perturbed 



where e > and — 1 < n < 0. Our model allows ex- 
act scaling during the radiation and matter dominated 
eras (provided by the first term) and drives the late-time 
acceleration in the recent epoch (provided by the second 
term) . Our model with the scaling initial conditions does 
not cause the numerical problem discussed in Sec. HI In 
the early radiation dominated epoch (e.g., starting from 
Oi/oo = 10^^^ in this paper), for small e = 10~^ ~ 10~^ 
the quantity ~ (1 + e)R'^ evolves with values slightly 
larger than unity so that Eq. ([T]) or ([2]) can be numeri- 
cally manageable within the precision of usual computing 
environment. The evolution of both background and per- 
turbed variables becomes numerically unstable as e gets 
smaller (e.g., e < 10~^ for Ui/ao = 10^^^). This is be- 
cause at the early epoch the quantity F is so close to 
the unity that the time-variation of F is not numerically 
manageable even in double precision accuracy. 

Although the detailed study is not given here, the value 
of e can be more tightly constrained by the solar system 
test. For example, according to the criterion given by 
Ref. [Hi, one expect e < IQ-^^ for R/H^ = 10^. As we 
explained earlier, such a small value of e is numerically 
problematic in the early era due to limited numerical 
precision currently available. Thus although the values 
of e we consider are likely to be too large considering 
the solar-system constraint, as shown below, the case of 
e = 10~^ and n = gives power spectra and perturbation 
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growth that are observationally indistinguishable from 
the predictions of the ACDM model in Einstein gravity 
(e = 0, n = 0). 

There is one numerical task needed during the back- 
ground evolution of f{R) gravity. Evolving Eq. ((1} or 
([2]) demands a fine-tuning process to satisfy the con- 
dition that the normalized Hubble parameter at the 
present epoch should be equal to unity, Hq = 1, where 
H = H/Hq. When the radiation and the matter energy 
densities at the present epoch arc fixed, one of the coef- 
ficients appearing in the f{R) functional form should be 
adjusted to match this condition for given initial condi- 
tions of R and i?', or vice versa. In our case we adjust 
q appearing in Eq. ((20)) to satisfy Hq = 1 for a given set 
of scaling initial conditions. The same numerical situa- 
tion also appears in the Einstein gravity with dark energy 
model based on a minimally coupled scalar field. 

In a modified gravity with the form f{R) = R + qR~^\ 
the method of imposing initial conditions for the back- 
ground evolution used in the literature is unclear to us. 
Furthermore, the coefficient q is not single-valued de- 
pending on the choice of n and initial conditions of R and 
R'. In an extreme case when n = —0.97, Ri — 10^^ Hq, 
and R'^ = at ai/ao = 10"""^^ with other parameters 
fixed with the Wilkinson Microwave Anisotropy Probe 
(WMAP) 7-year best-fit parameters (see below), there 
are nine multiple values of q satisfying Hq = 1 in the 
range -0.7222 < q/HQ^^"^ < 0. Under the same (non- 
scaling) condition, multiple q values are also obtained in 
our /(i?) gravity. On the other hand, our initial con- 
ditions for the scaling evolution presented in Sec. Illll are 
exact and general so that they can be applied to any dom- 
inant fluid with a constant equation of state w. With the 
scaling initial conditions, q is single- valued over the whole 
range of n considered. 

In our analysis we considered cases where the behavior 
of our f{R) gravity, Eq. (|20p can be approximated by the 
first term (i?^"*"*^; higher power of R) at early times and 
is dominated by the second term {qR~", — 1 < n < 0; 
lower power of R) at late time during the acceleration 
phase. However, it should be emphasized that the as- 
sumption that the higher power of R should dominate 
at early epoch and the lower power of R should domi- 
nate at later epoch is not necessarily the case due to the 
nonlinearities of the f{R) gravity theory [l^]. Because 
the scaling background evolution guaranteed by the scal- 
ing initial condition leads to decreasing R value in time 
for u; > — 1 [this follows from Eq. (|14p ]. this allows the 
transition from higher to lower i?-power regimes in our 
case. 

Figure [1] shows background evolution of our f{R) grav- 
ity with e = 0.001 and varying n (color solid curves). 
The background parameters shown are density parame- 
ters fii, energy densities /Xj, equation of state wx, rela- 
tive Hubble parameter (H/Hacdm) and distance mod- 
ulus {p — /xacdm) with respect to the fiducial ACDM 
model. As a fiducial model we take a flat ACDM uni- 
verse with parameters Hmo = 0.274 (ilco = 0.2284 and 
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FIG. 1: Background evolution in the f{R) = R^+' + qR'"^ 
gravity with e — 0.001 and varying n = —0.01, —0.05, —0.10, 
-0.15, -0.20, -0.25, and -0.30 (color solid curves). Results 
for e — 10~^ and n — —0.1 are added with blue dashed curves. 
(Top panels) Evolution of fii and Hi as a function of scale 
factor a{t), where i = R, M, X indicates radiation, matter 
(baryon plus CDM), and gravity sector, respectively. (Mid- 
dle and bottom panels). Evolution of Q,x, equation of state 
of X-component wx, Hubble parameter H / Hacdm, and dis- 
tance modulus A/i — /i — /iACDM relative to the ACDM model, 
where p. represents the distance modulus. In all panels ACDM 
predictions are shown as thick black curves. The Hubble pa- 
rameter data is taken from Ref. [l^, and the binned SNIa 
data points are from the UNI0N2 sample [23| . 



^ 0.0456), Qao = 0.7278, h = 0.704, = 0.963, 
erg = 0.809, To = 2.725 K, = 0.24, = 3.04 with 
reionization optical depth t = 0.087 based on the WMAP 
7- year observations [13]. In all f{R) gravity models, we 
set Qxo = ^^Ao- With the fiducial model parameters, the 
e has an upper bound ebbn = 0.011314 so that the ini- 
tial contribution from the dark energy is lower than the 
maximum amount allowed by the big bang nucleosynthe- 
sis (BBN) calculation, Qxi < 0.045 [2l|. We also added 
blue dashed curves in Fig. [T] to represent the background 
evolution for e = 10^^ and n = —0.1, demonstrating that 
the energy density of AT-component (here Ffix) becomes 
less dominant than the dominant fluids {^n and /xm) as 
e gets smaller. The background observables (Hubble pa- 
rameter and distance modulus) for e < lO"'^ are largely 
insensitive to e and are very similar to those for e = 0.001. 

In order to calculate the matter and CMB power spec- 
tra, we solve a system composed of matter (baryon [b] and 
CDM [c]), radiation (photon [7] and neutrinos [i/]), and 
the f{R) gravity sector as dark energy. The perturbation 
in the radiation component has been handled by using 
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FIG. 2: Top panels: baryonic matter (left) and CMB tem- 
perature anisotropy (right) power spectra in the f{R) — 
R^^" + qR^" gravity models with n = —10^^ and varying 
e = 10"^ 10"^ 5xl0"^ 10"^ 5xl0"^ lO^"*, and 10"^ Bot- 
tom panels: The same as in the top-panels but with e = 10~^ 
and varying n = -10'^, -10"*, -10~^ -0.01, -0.1, -0.3, 
and —0.5. The results for e — 10^^ and n = are shown as 
brown dashed curves, which are very similar to the ACDM 
predictions (black curves). In all models including ACDM 
model, the same initial amplitude has been assumed. The 
ratios of f{R) gravity power spectrum to ACDM prediction 
are also shown in the bottom region of each panel. For mat- 
ter and CMB power spectra, recent measurement from the 
Sloan Digital Sky Survey (SDSS) DR7 luminous red galaxies 
(LRG) [3 and the WMAP 7-year data (including the cos- 
mic variance) [2^ have been added (grey dots with error bars) 
together with fractional errors of observed spectra. 



the Boltzmann equation or tight coupling approximation 
(see Ref. [1^ for our set of equations and the numerical 
methods). As the initial conditions of background and 
perturbation variables we use Eqs. ([T ^ - p^ and 
with w = ^ and f2„,o = ^ro- All the perturbed variables 
we use are spatially gauge invariant [2j]. We solved the 
system by adopting the CCG as the temporal gauge con- 
dition. The CCG is the same as the synchronous gauge 
without the gauge mode. For the matter power spectrum 
we present the power spectrum of baryonic density per- 
turbation based on the CCG which is a gauge-invariant 
concept. The CMB anisotropy is naturally gauge invari- 
ant. Figure [5] shows the matter and CMB temperature 
anisotropy power spectra for n = —10~^ with varying e 
(top) and for e ~ 10~^ with varying n (bottom panels). 
We omit CMB polarization power spectra. Note that the 
baryonic matter power spectra is very sensitive to both 
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FIG. 3: Top-left: evolution of density perturbations Si (i = 
7, 6, c, X) at comoving scales k = 10~^ and 10~^ hMpc~^ in 
f{R) = R^+' + qR~" gravity with n = -10"^ for e = 0.001 
(solid) and e = 10^^ (dashed curves). Perturbation growth 
of photon (7), baryon (6), CDM (c), and 5R/R are repre- 
sented as red, yellow, blue, green curves, respectively. Other 
panels: evolution of growth factor g = 5b /a (normalized to 
unity at present) at comoving scales k = 10"^, 10~^, and 
10^^ /iMpc^^ for models with parameters used in the top- 
panels of Fig. [2] (n = — 10~^ and varying e) with the same 
colored code. Growth factors of the ACDM model are shown 
as black curves. Note that since we consider interactions be- 
tween radiation and baryon components, the ACDM growth 
factor is mildly scale-dependent. Yellow dots with error bars 
indicate the ACDM growth factor expected from the future 
X-ray and weak-lensing observations [2a ]. 



e and n while the CMB power spectra (including polar- 
ization) are mildly sensitive to the variation of n and are 
insensitive to e. 

Our results for e = 10^^ and ?i = (shown as brown 
dashed curves in the bottom panels of Fig. [2]) are very 
similar to the ACDM predictions (black curves); the de- 
viations are less than 2% level at scales k < 0.2 /iMpc~^ 
for matter power spectrum and less than 0.5% level for 
CMB power spectrum. Our calculations of power spec- 
tra for e = 10^^ and nonzero n are consistent with those 
presented in Ref. [6|] that considered f[R) = R + qR~"' 
gravity. We have also checked that the power spectra for 
e = 10^'' are indistinguishable within the 1% precision 
from those obtained by evolving the perturbed variables 
under the model which uses the ACDM model in Einstein 
gravity at the early epoch (a/ao < 0.001) and suddenly 
switches to the /(i?) = R + qR~" gravity thereafter. 

Figurc[3]displays the perturbation growth of individual 
components {S^, Sb, Sc, 6R/R) and baryon perturbation 
growth factor g = Sb/a at recent epoch in our f{R) grav- 



log a(l) 



_ 1 1 . 1 1 . 1 1 . 1 1 1 1 1 _ 
- k-10-2 h Mpc-' - 

^ = 






" 






FIG. 4: Same as Fig.[3]but for /(i?) gravity models (e = 10"^ 
and varying n) used in the bottom-panels of Fig. [5] with the 
same colored code. 



ity with n ~ — 10~^ and varying e. We see that each 
perturbed variable foUovifS with each other at the early 
epoch, which demonstrates the scaling evolution. In the 
growth factor panels we have added ACDM mock growth 
factor data points that are expected from the future X- 
ray and weak-lensing observations [IBl- The mock data 
points, located in the redshift range z = 0-2 with equal 
interval of Az = 0.2, have been generated by adding 1% 
random noise to the ACDM growth factor at comoving 
scales k = 0.001, 0.01, and 0.1 /iMpc"^ The perturba- 
tion growth factor at the small scale (fc = 0.1 /iMpc~^) 
is very sensitive to the variation of e. It is also sensitive 
to the variation of n, which is shown in Fig. 21 

From Figs. [2]|4] we notice that the CMB power spec- 
trum is less sensitive to the model parameters than the 
baryon density power spectrum and the baryon growth 
rate. Thus, in the following we use the baryon density 
power spectrum and the baryon growth rate together 
with the SNIa data to constrain the model parameters. 
Figure [S] shows likelihood distributions of /(i?) gravity 
parameters constrained by the recent UNI0N2 SNIa data 
set (without systematic errors) [lOl, the SDSS DR7 LRG 
power spectrum (PS) [IJl, and the ACDM mock growth 
factor data at the comoving scale fc = 0.1 /iMpc^^ (in 
the bottom-right panels of Figs. |3] and |4]). Our like- 
lihood estimation is tentative since we have explored 
only the e-n space while other cosmological parame- 
ters are fixed with the WMAP 7-year best-fit parame- 
ters. The likelihood distribution has been calculated from 
L — Lniax cxp(— Ax^/2), where Lmax has been taken as 
the maximum value within the parameter space we have 
explored, 10"^ < e < ebbn and -0.5 < n < -10"^. 




FIG. 5: Likelihood of f{R) gravity parameters constrained 
by the UNI0N2 SNIa sample (blue contours), the SDSS DR7 
LRG power spectrum (red contours), and the ACDM mock 
growth factor data at fc = 0.1 /iMpc~^ (green contours). The 
levels of contours have been determined by setting Ax^ ~ 2.30 
(solid), 6.17 (short-dashed), 11.8 (long-dashed curves) relative 
to the x^-minimum, mimicking la, 2a, 3a confidence levels, 
respectively. Arrows indicate directions to which the likeli- 
hood increases. Triangles and squares indicate the models 
shown in top and bottom panels of Fig. [2] (Fig. [3] and Fig. [4] 
for the perturbation growth rate), respectively, with the same 
colored code. 



The 1(7, 2(7, 3(7 confidence levels have been roughly de- 
termined by setting Ax^ = 2.30, 6.17, 11.8 from the x^- 
minimum, which is valid only for the Gaussian distribu- 
tion. In the x^-estimation with the SDSS LRG power 
spectrum, we consider the convolution effect caused by 
band-power window functions and exclude data points 
where the non-linear clustering dominates (see Appendix 
C of Ref. [131 for detailed description). As shown in Fig. 
m f{R) gravity parameters, e and n, are very sensitive to 
the growth factor at small scales, and are already tightly 
constrained by the current measurement of galaxy power 
spectrum. For ACDM mock growth factor data, the 
likelihood distribution suggests that only the parameter 
space extremely close to the ACDM model is allowed. 



V. DISCUSSION 

In this paper we have studied a /(i?)-gravity based 
dark energy model with early scaling era. We have pre- 
sented initial conditions of background and perturbed 
variables during the early scaling evolution regime in the 
modified gravity with a pure power-law form f{R) = 
in the early era. With these initial conditions, the 
modified gravity with a form f{R) = R^'^'^ + /de(-R) 
where the second term drives the late-time acceleration 
becomes free from the numerical difficulty that is usu- 
ally confronted during the background evolution in the 
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early radiation dominated era for f{R) = R + /de(^) 
gravity. Our initial conditions are general so that the 
scaling density evolution of the X-component is assured 
for any dominant fluid with a constant equation of state 
parameter w. 

As a possible dark energy model we have considered 
the gravity with a form f{R) = R^'^'^ + qR~^ and com- 
pared the evolution of the background and perturbation 
variables in this gravity with the recent observational 
data and the ACDM mock data. The present observa- 
tional data already severely constrain our model param- 
eters n and e so that only parameters extremely close 
to the ACDM model is allowed. We found that the 
power spectrum of baryon component and the pertur- 
bation growth factor at small scales are more sensitive to 
the f{R) gravity parameters than the SNIa distance mod- 



ulus and the CMB anisotropy power spectra (see Figs. [Il- 
ls]). Therefore, precise measurement of the perturbation 
growth is essential to tightly constrain our f{R) gravity. 
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